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Abstract 

, We study the relationship between the holomorphic unitary connection of Chern- 

\^ ' Simons theory with temporal Wilson lines and the Richardson's exact solution of the 

^SJ , reduced BCS Hamiltonian. We derive the integrals of motion of the BCS model, their 

' eigenvalues and eigenvectors as a limiting case of the Chern-Simons theory. 

o 



Oh! 

r-| ! 1 Introduction 

> : 
•1—1 , 

^ ' Chern-Simons theory has emerged in the last years as a very useful field theoretical 

' model for the description of interesting phenomena in Condensed Matter Physics. It plays 

a leading role in the description of the quantum Hall effect . It has been also advocated 
to be connected with hight Tc superconductivity Q. In this paper we will show that it 
also appears to be closely related to the Bardeen, Cooper and Schrieffer (BCS) theory of 
superconductivity Q. 

One of the corner stones in Condensed Matter Physics is the pairing model of Bardeen, 
Cooper and Schrieffer (BCS) which explains the properties of "low Tc" superconductors Q 
and several nuclei The BCS model was solved in the Grand Canonical Ensemble by a 
variational ansatz, which is asymptotically exact in the number of fermions forming pairs. 
A reduced version of the BCS model, where the pairing interaction has the same strength 
between all pairs, was solved exactly in 1963 by Richardson in the Canonical Ensemble, 
where the number of Cooper pairs is fixed fS], ^, ^. 

Despite of the interest of this exact solution it escaped the attention of the condensed 
matter and nuclear physics communities, until recently with the advent of ultra-small 
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superconducting grains |^] , whose study is based on the BCS model with a fixed number 
of pairs (for a review see ) . 

Partially motivated by Richardson's work, in 1976 Gaudin proposed a family of com- 
muting spin Hamiltonians whose exact diagonalization paralleled Richardson's solution 
[10|. However a relation between the two models was not clear at that time, despite of the 
fact that in the limit where the BCS coupling constant g goes to infinity, the Richardson's 
solution turns into the Gaudin's one. 

A first step to understand the relation between the Gaudin's spin Hamiltonians and the 
reduced BCS Hamiltonian came with the work of Cambiaggio, Rivas and Saraceno (CRS) 
who proved the integrability of the reduced BCS model in terms of a family of commuting 
operators, which are nothing but the Gaudin's Hamiltonians plus a perturbation propor- 
tional to the inverse of the BCS coupling constant g |11]. CRS also constructed the BCS 
Hamiltonian as a linear combination of the perturbed Gaudin's operators. Unfortunately 
these authors were unaware of Richardson's and Gaudin's works and consequently they 
did not obtained the eigenvalues of the BCS conserved quantities. 

The eigenvalues of the CRS operators were given in the reference p!^ , using Conformal 
Field Theory (CFT) methods, and generalized to other models in |l^ The aim 
of reference ||l^ was to place the Richardson's exact solution and its integrability in the 
framework of CFT, relating it to the Wess-Zumino-Witten models with Kac-Moody alge- 
bra SU{2)k- Using CFT methods it was shown that the Richardson's wave function are 
obtained as conformal blocks of a perturbed S'J7(2)fc-WZW model when the level k goes 
to a critical value —2 |12]. 

The CFT interpretation of the BCS model turned out to be closely related to the 
work of Babujian [16|, who in 1993, used the so called off-shell algebraic Bethe ansatz, 
to re-derived Gaudin's exact solution |]l^, |l^]. Babujian also observed that the Gaudin's 
eigenstates can be used to build the conformal blocks of the WZW models. 

The similarity between the works |jl^ and [16, ^] suggested that the Richardson's 
solution of the reduced BCS model should also be derivable using the off-shell Bethe ansatz 
method. This was done by Amico, Falci and Fazio pi] and later on clarified in references 



[19, 20, 21 [, where the BCS coupling constant parametrizes a boundary operator that 
appears in the transfer matrix of the inhomogenous vertex model, whose semi-classical 
limit gives rise to the CRS conserved quantities. 

The old results by Richardson and Gaudin's, concerning the norm of the eigenstates 
and the occupation numbers [22, |2^, has been generalized to other operators in references 
[24, 1£[ using the "determinantal" techniques developed in [25, [2^ . 

The outcome of all these works has been to clarify the integrability and relationship 
between the reduced BCS and Gaudin's model at least at a formal level. There remains 
however the question concerning the "geometrical" or field theoretical origin of the integra- 
bility of BCS. For the Gaudin's Hamiltonians this is given by the underlying WZW model, 
and ultimately by the Chern-Simons theory [^]. The BCS model is related to some sort of 
chiral perturbation of the WZW model, characterized by modified Kniznik-Zamolodchikov 
equations [||, |l|, ||. 

In this paper we shall show that the field theoretical origin of BCS can be traced back 
to a SU{2) Chern-Simons (CS) theory interaction with a one-dimensional distribution 
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of coloured matter which breaks both gauge and conformal invariance. This connection is 
quite remarkable because Chern-Simons theory has been advocated to be mainly connected 
with effective descriptions of fractional quantum Hall effect and high superconductivity, 
but never with standard BCS superconductivity. The Chern-Simons theory is not defined 
in the physical space, which might be two or three-dimensional, but rather in the complex 
energy plane which is always two-dimensional. This explains why this field theoretical 
connection of BCS theory remained unveiled for so long time. The connection of Chern- 
Simons theory with BCS model can be understood in a more general framework when we 
consider a scaling limit of the twisted Chern-Simons theory defined on a torus, that is, the 



twisted elliptic Chern-Simons theory. On a torus the KZ equations [28| are replaced by the 
Knizhnik-Zamolodchikov-Bernard (KZB) equations which depend on the coordinates 
Zn of the punctures, the moduli of the torus r and a set of parameters uj characterizing 
the toroidal flat gauge connections. The later parameters uj define the twisted boundary 
conditions for the WZW fields on the torus. 

The main result of this paper is to show that, for a generic simple simply connected, 
compact Lie group G, the Richardson equations, the CRS conserved quantities and their 
eigenvalues arise from the KZB connection and their associated horizontal sections. This 
will be done in a limit where the torus degenerates into the cylinder and then into the 
complex plane. In this limiting procedure the generalized BCS coupling constants appear 
as conjugate variables of the parameters uj, when this parameters go to infinity. This gives 
the G-based BCS models a suggestive geometrical and group theoretical meaning. 

The organization of the paper is as follows. In section 2 we briefly review the reduced 
BCS model and its exact solution. In section 3 we consider a perturbation of the CS model 
defined on the plane, we derive the KZ equations and study their connection with the CRS 
conserved quantities and the BCS model. In section 4 we introduce the KZB connection 
related to the CS theory on the torus. In section 5 we derive the Bethe ansatz of the 
integrable Hamiltonians introduced in sections 3 and 4. Finally we state our conclusions. 
We have also included an appendix with particular examples of the general equations 
contained in the main text. 

2 Review of the exact solution of the BCS model 

The reduced BCS model is defined by the Hamiltonian 12 1 



Hbcs= Snaci^c^^ - gdJ2ci+cl_c^>_c^>+ (1) 

n,a=± n,n' 

where c„^± (resp. cjj_|_) is an electron destruction (resp. creation) operator in the time- 
reversed states \n,ziz) with energies e„, d is the mean level spacing and g is the BCS 
dimensionless coupling constant. The sums in (|l]) run over a set of doubly degenerate 
energy levels en(n = 1,... ,A^). We shall assume that the energy levels are all distinct, 
i.e. Em 7^ En for m ^ n. The Hamiltonian (||) is a simplified version of the reduced 
BCS Hamiltonian where all couplings have been set equal to a single one, namely g. This 
is the model that is commonly used to describe ultrasmall grains and accounts for the 



3 



scattering of pairs of electrons between discrete energy levels that come in time-reversed 
states. Hereafter we shall refer to (|l]) simply as the BCS Hamiltonian. 

Richardson had long ago solved this model exactly for an arbitrary set of levels, e„ 
1^, 0. To simplify matters, we shall assume that there are not singly occupied electronic 
levels. As can be seen from these levels decouple from the rest of the system; they are 
said to be blocked, contributing only with their energy £^ to the total energy . The above 
simplification implies that every energy level n is either empty (i.e. |vac)), or occupied 
by a pair of electrons (i.e. cj^ _,_cjj _|vac)). Denote the total number of electrons pairs by 
M. Then of course M < N . The most studied case in the literature corresponds to the 
half-filled situation, where the number of electrons, N(, = 2M, is equal to the number of 
levels |9|. In the absence of interaction (i.e. g = 0), all the pairs occupy the lowest 
energy levels forming a Fermi sea. The pairing interaction promotes the pairs to higher 
energies and eventually, for large values of M, all the levels are pair correlated, giving rise 
to superconductivity 

In order to describe Richardson's solution one defines the hard-core boson operators 

bn — Cn,—Cn,-\-, = c|^^_|_cjj _, Nn = b^bn (2) 

which satisfy the commutation relations, 

K,bl]=5n,n' il-2N^) (3) 
The Hamiltonian (||) can then be written as 

HbCS = J2'^^ribib^-gJ2 ^n^n' > (4) 
n n,n' 

where we have set d = 1 (i.e. all the energies are measured in units of d). Richardson 
showed that the eigenstates of this Hamiltonian with M pairs have the (unnormalized) 
product form js], ^, ^ 

M N I 

\M)r = X{B,\w^c), B, = Y, _ bi (5) 

j=l n=l 

where the parameters Ej (j = 1,...,M) are, in general, complex solutions of the M 
coupled algebraic equations 

^ .^ Ef - ~ ^ 2en - Ei ' 

^ 3'=H¥^]) ^ J n=l " J 

which are a sort of Bethe ansatz equations for this problem. The energy of these states is 
given by the sum of the auxiliary parameters Ej, i.e. 

M 

£{M) = Y,E, (7) 
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The ground state of Hbcs is given by the solution of eqs.(^) which gives the lowest value 
of 8{M). The (normalized) states (^) can also be written as [^] 

^ N 



where the sum excludes double occupancy of pair states (rij ^ nj for i ^ j) and the wave 
function if; takes the form 



M 1 

V'^K,...,n,,)= y: u^, — ^ (9) 

VeSM k=i ^^"^ - ^^fc 



The sum in runs over the set of permutations Sm, of 1, • • • , M. The constant C in (^) 
guarantees the normalization of the state |Q (i.e. r{N\N)ji = 1). 

A well known fact about the BCS Hamiltonian is that it is equivalent to that of a XY 
model with long range couplings and a "position dependent" magnetic field proportional 
to En- To see this let us represent the hard-core boson operators (^) in terms of the Pauli 
matrices as follows, 

6„ = a+, 6], = cr^, Af„ = ^(1 - 0-2) (10) 

in which case the Hamiltonian ^ becomes 

N 

Hbcs = HxY + Y.^n + g{N/2 - M) 

N (11) 

HxY = -Y. 2enr° - §(r+ T- + T- T+) 

n=l 

where the matrices 

N 

T'' = Y.< (« = 0,+,-) 

n=l (12) 

^0 ^ I z + ^ + - ^ - 
' n 2 rtJ 'n n 1 'n 

satisfy the SU (2) algebra, 

[T'',T^] = ffT" 
/+o = /o- = _i, /+- = 2 
whose Casimir is given by 

T-T = r°r° + -(T+r~ + r~r+) (i4) 

2 



(13) 
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This spin representation of the BCS model is the appropriate one to study its inte- 
grabihty, as was shown by Cambiaggio, Rivas and Saraceno by exphcitly constructing a 
set of operators, Rn {n = 1, . . . , A^), which commute with H^cs &iid moreover allow a 
reconstruction of the BCS Hamiltonian which is given by Hbcs = gJ2n'^^n Rn + ctes. 
In the next sections we shall generalized this construction to arbitrary Lie groups and 
representations. 



3 Chern-Simons Theory and BCS Hamiltonian. 

Let G be a simple, simply connected, compact Lie group. The Chern-Simons theory of 
N interacting heavy coloured particles (temporal Wilson line insertions) located at points 
zi, Z2, ■ ■ ■ , zn of the 2D plane C is defined by the action pO] 



S{A,gi,...,gN) = ^ J ti(^AAdA + ^AAAAA^ 

N roo ^ (15) 

+ dt {An,gn\t)Ao{zn,t)gn{t)+g-\t)-gn{t)) 

where j4 is a G-gauge field, gn a G-valued chiral field describing the dynamics of the 
coloured particle located at the insertion point z„ and A„ the highest weight of its G- 
colour representation. (• , •) denotes the Killing scalar product in the Lie algebra g of G 
normalized so that long roots have length 2, we use this scalar product to identify Q and its 
dual g*. The quantization of the theory (^) yields quantum states which have to satisfy 
two constraints. In the Bargmann representation they are holomorphic functionals ^{A^) 
of the (0, 1) component of the gauge field A^ = ^{Ai + 2^2) and take values in the tensor 
product of the representation spaces <^Va^ of the insertions. The remaining constraint is 
Gauss' law 

^^-'tWT^^^--) = --dzAl^iA,) + tl)6{z - Zn)^{A,), (16) 



M|(z) 



n=l 



where t"" is an orthonormal basis of the Lie algebra g of G and t^^-^ denotes the operator 
representation on of t°- S g. The law (|T^ ) governs the dependence of the quantum 
states ^ under a G'^-valued complex gauge transformation g, i.e. 

N 

^{g-'A,g + g-^d,g) = e^Swzw(A,,9) J| g(^„)-i ,f (^,) (17) 

n=l 



where 



Swzw{A2,g) = J dzdztvg ^{dzg)g d ^ti{g ^dgf 

+ J dzdz tiA^gdzg'^ 



is the Wess-Zumino-Witten action. In equation (|17|) denotes the operator representa- 
tion on Va^ of the element g of G*^ . The global implementation of this dependence requires 
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that the couphng constant k is an integer. Since any gauge field A can be rewritten as 
A = g~^dzg in terms of a complex gauge transformation g, all physical states are deter- 
mined by its value at the trivial gauge field configuration A^ = 0, i.e. ^'(0). On this 
configuration, the infinitesimal dependence of the physical states on the gauge field can be 
obtained p^] from Gauss' law (p!6|) 



vr 



6A1iz) 



nA.)i__, = j:-^m (18) 

n=l " 



and 



5AUz)5AUw) ^''a=o {z-wf ^ >^^^iuj-z^){z-z^) 

,^^{W-Zn){z-Zm) 

where the sum in the repeated index a is implicitly assumed and l/{z — w) is the Green 
function of the operator ^dg- From these expressions and the Sugawara form of the stress 
tensor 



7r2 , r S"^ kdimG 



2{k + h"^) ' 5Al{z)6Al{z) ' 2{k + h'^)w^zl5Al{z)6Al{w) {z - w)^ 



(20) 



it is straightforward to show that the dependence on the translation of the insertion at z„ 
is governed by the Knizhnik-Zamolodchikov equation [28] 

5.„*(0) + i E ^^^^^(0) = 0, (21) 



where k = k + and hy is the dual Coxeter number. In the same way the dilation 
transformation of the plane implies that 



N ^ 



z„9,„^(0) + -Hg^{Q) = (22) 



1 ^ 

n=l 

where the operator 

in the second term of the equation is essentially the coupling term in the BCS Hamiltonian 
(eq ([TT|)). This suggest that a possible slight modification of the Chern-Simons theory can 
trace back the geometric origin of the BCS Hamiltonian. 



^ The metric dependence of the theory is induced both by the Bargmann normalization of quantum 
states (k) and the removal of spurious gauge degrees of freedom (/i"^) 



7 



Let us consider the perturbation of Chern-Simons theory by a one-dimensional distri- 
bution of coloured charges. The effect of this perturbation is encoded by a new interaction 
term of the form 



where ^ is any element of the complexified Lie algebra g , C is a non-self-intersecting closed 
curve of C which encloses all point-like insertions Zn of heavy coloured particles and where a 
one-dimensional .^-coloured charge distribution is localized. Because of its ^ dependence this 
perturbation is not gauge invariant under gauge transformations not completely included in 
the normalizer of ^. On the other hand, conformal invariance is broken by the special choice 
of the closed curve C. However, the full Chern-Simons theory defined by S{A,g) + S^{A) 
can be quantized in a similar way and the corresponding quantum states are holomorphic 
functionals ^[A^) taking values in the same vector space ®Vf^^ because the new insertions 
have one-dimensional representations. The essential difference is encoded by the Gauss' 
law which now reads 






(24) 



This implies that the variation (|T9|) now becomes 




where ^(n) denotes the action of ^ on the representation A„ of the coloured particle sitting 
at Zn- Using the Sugawara form of the stress tensor (|^) it is easy to show that the physical 
states must satisfy a modified Knizhnik-Zamolodchikov equation 




(26) 



In a similar way the dilation transformation now yields the equation 




(27) 



and, what is remarkable, in the case G = SU{2) the operator 




(28) 
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does coincide up to a factor 1/ g and additive constant with the reduced BCS Hamiltonian 

The solvabihty of the reduced BCS Hamiltonian worked out by Richardson can now be 
understood in this geometrical approach from the existence of a complete set of operators 

Rn=Y.'^+kn) (29) 

commuting with = 0, as first pointed out in Ref. [11|. The operators i?„ also 

commute each other = 0) and provide a complete set of constants of motion. 

This remarkable feature generalizes the well known integrability of the Hamiltonian (p^) 
which commutes with the Gaudin Hamiltonians. The only difference is due to the presence 
of the 1^ term in (^). The commutativity of the Hamiltonians i?„ also means that the 
connection (Knizhnik-Zamoldchikov connection) defined by them in the bundle of Chern- 
Simons states over the space of insertions R-^ is flat. In fact, this flatness property of 
the KZ connection is not only verified on the bundle of CS states but in a more general 
bundle of states where the extra condition of global colour invariance imposed by Gauss' 
law is not required. This colour neutrality condition comes from the integration of Gauss' 
law equation (|^) and establishes that CS states can only take values on the subspace 
C (^V\„ of invariants under the diagonal action of the Cartan subgroup of G (here and 
below we take ^ in the Cartan subalgebra of G) and as we shall see later it will correspond 
to the half-filling regime of BCS superconductivity. In fact, this restriction can be ignored 
if we include an extra charge insertion zq at very large distance of all the other insertions 
with the only purpose of neutralizing their colour charge. In this way, the whole BCS 
theory emerges from CS in the singular limit k ^ as we shall see in next section. This 
relation existing between Chern-Simons theory and the BCS theory can be even better 
understood from a elliptic point of view. 

4 Knizhnik-Zamoldchikov-Bernard connection. 

Let us consider now the Chern-Simons theory on a compact two dimensional space 
with genus one. This theory was studied in [^, Let us summarize the main results. 
To quantize the theory one has to pick up a complex structure in the torus C/ (Z + rZ) 
with T2 = Imr > and a local systems of coordinates Zn for the n*^ Wilson line insertions. 
Then, quantization proceeds in a similar way as for the plane or the sphere. In Bargmann's 
quantization the states are holomorphic functionals of the (0, 1) part of the gauge 1-form 
with values in the tensor product of the representation spaces ®V\^. The states are subject 
to the Gauss' law that once integrated tells us how the states transform under complex 
gauge transformation. Every (generic) gauge field can then be gauged out to a flat one of 
the form Au = 7rudz/T2 with u G f)*^, the (complexified) Cartan subalgebra of G, and the 
states are effectively described in terms multidimensional theta functions 6{u) with values 
in the subspace C ®Va„- The states are also subject to extra regularity conditions. See 



H, |3|, m for details. 

The dependence of the quantum states on the positions of the point-like insertions can 
be derived in a similar manner as in the planar case. However, in the elliptic case Gauss' 
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law ( |T6|) does not provides in a simple manner the infinitesimal variation of quantum states 
with respect to arbitrary gauge fields. In fact, the required Green function of the operator 
Dz = dz + t^u/t2 does depend on u. The trick is to remove the u dependence from the 
differential operator Dz = dz + t^u/t2 by means of a singular gauge transformation 



N 



n=l 



The u dependence is traded in this manner by a non-trivial twisted boundary condition 
for Chern-Simons states and the Green functions of the operator dz- We further extract a 
normalization factor accounting for the volume of the complex gauge orbits 



rankG . , / 2niTl+2ni(u,a) 



o>0 1=1 a 

with the product in a (resp. a > 0) extending to the set of roots (resp. positive roots) of 
G. If we then define 9 by 

7(tx) = Ii{u)e{u), 

the dependence of the CS quantum states on the coordinates of the point-like charged inser- 
tions and on the elliptic modular parameter r |35] can be expressed in a simple geometrical 
meaningful form: Chern-Simons states describe horizontal sections 

V^6l(n) = 

with respect to the flat Knizhnik-Zamolodchikov-Bernard (KZB) connection A defined by 
=0, A-z^=Q, Ar=\Ho{T,z), Az„= ^Hn{T,z). (30) 

with 

N , 

H ( -2^5^(T(„,a)(^n-2m) (ea)(„)(e-a)(^) 

(31) 



47r Svr 



n,m=l 



j=l m^n ^ " 



(32) 



on the vector bundle of Chern-Simons states over the moduli of punctured elliptic curves. 
In equations (|3l|)(p^) p and ax denote the Green functions of ^dz that can be expressed 
in terms of the Jacobi theta function 

^ '<^l^_-^y^7Ti(l + l/2)^T+27Tiz(l + l/2) 
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by 

and 



'd{x)'d{z) 



and we use the notation of Refs. Cq, represents the step generator associated 

to root a of g and {hi,i = 1, . . . , r} is an orthonormal basis of the Cartan subalgebra, so 
that u = J2j '^jhj- 

A natural hermitian structure in this bundle was introduced in Ref. [^9| by means of 
the scalar product a la Bargmann, in that paper it was also shown that KZB connection 
is the unique (1,0) unitary connection in this bundle. 

It is easy to check that the KZB connection leave invariant and restricted to is 
flat for any value of k. This implies, in particular, that Hamiltonians Hn restricted to 
commute for n = 0, 1, . . . , and determine then a quantum integrable system that admits 
a family of eigenvectors using the appropriate Bethe ansatz [^]. Below we shall discuss 
the solutions to the problem of eigenvalues. 

Notice that in this case the equation corresponding to a global dilation is more involved 
that in planar case which makes the connection with a Hamiltonian of BCS type less 
transparent. However, we shall show that in the limit r — > ioo where the torus degenerates 
into the cylinder we shall recover a connection with the BCS theory. In particular, after 
rescaling t rjr and taking the limit — > oo we obtain the new family of commuting 
Hamiltonians, Hq and Hn that read: 

47r 4 ^sm (7r(u, a)) 



i=l mT^n (-34-) 

+ 7r 51 Z!c°*(^('"'"))(^")(n)(e-a)(,„) 

with e„ = En(ea)(„)- 

Next we take the limit from the cylinder to the plane by sending variables z„ and u to 
the origin, we rescale Zn 1— > (zn and u 1-^ (u, Rq = C'^Hq, Rn = C-^n- In the limit of ^ — > 
the new commuting Hamiltonians are: 



Mz) = 2^A„ - E , (35,a) 



47r " A-K ^ (u, aV' 
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We finally look for solutions to the eigenvalue problem for the Hamiltonians (|35| ) of the 
form 

T(n,i) =e<«'">F(u,z) (36) 

such that ^ G f)*^ and F{u/e,z_) is for generic u and z holomorphic in e G C, i. e. we assume 

F{u/e,z) = Fo{z)+eFi{u,z) + ... (37) 

If solutions of this form exist (later we shall see that this is the case) then Fq is an 
eigenvector of Hamiltonians, Rn = lim„_+oo e~^^'"^i?rae^^'"\ and 



In the case of SU{2) these are the conserved quantities build by CRS The re- 

duced BCS Hamiltonian is essentially given hj H = J2n ^nRn (^) and hence it obviously 
commutes with all Rn- 

5 Bethe ansatz and horizontal sections. 

We shall address now the question of finding solutions to the equations that appear in the 
different limits. Following l3^ we define for /? G f)* 



UJ/3{z) = (Ti^u^p){z)dz + —dx(7(^u,l3){z)dT 

note that a; is a closed form in z, r space. 

Now decompose the sum of highest weights in simple roots (we assume that the de- 
composition exists): 

K N 
s=l n=l 

and for any partition I£ = {Ki, • • • , Kn) of K in N non-negative integers, relabel the roots 
a = (ai, . . . , ax) = ("i,i, ■ ■ • , «2,i; • • • , «Af,_ft:jv)- Introduce a complex variable ys for 

any simple root as and denote y = {yi,..., yx) = (yi,i, • • • , yi,Ki,y2,i, UN^Kf^)- 
Set the form 

^K,a{T,u,z,y) = a;^i_,(yi,i - 2/1,2) A 0^/3, 3(2/1,2 - ^1,3) A • • • A W/j^^^^ (yi^/^i - zi) 
A 

(39) 

A ^Pn,! iVN,! - yN,2) A u;/3jv,2 {yN,2 - yN,3) A • • • A {yN,KN - Zn) 

®n (e_„„^ •••e_a„,^^)(„)|A) 

where (3n,i = Yl\'=i '^n,i' and |A) = (8)„|A„) stands for the tensor product of highest weight 
states of representations A^. 
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And finally define 

^{T,u,z,y) = Y^ {-l)\''^\^K,VgSj,u,z,Vy) 
K V&Sk 

The form satisfies the equation 

95 A 17 + dr A H^n + Y,dzn^ H„Q = (40) 



where S is the multivalued holomorphic function 

(41) 



S{T,z,y) = ^ (A„, A„/) ln7?i(z„ - z„/) - ^(A„, a^) lni?i(z„ - 

n<n' _ n,s 

+ Y{as,as')ln'di{ys - Vs') 

s<s' 



with ^i{z) = 'di{z)/i!)'i{0), and d stands for the holomorphic differential in t,z and y 
variables. 



Equation (40) was stated in |38] and plays the key role in order to get horizontal sections 
of the KZB connection (|30| ) from integrals of form in cycles with coefficients in a local 



la 

system determined by e^'^. See ||38| for more details. 

We are interested in a complementary aspect of eq. ( |40| ) that relates it directly to the 
Bethe ansatz solutions of Hamiltonians Hi. Let the function G{t,u, z,y) be given by 

G dr A d^z A d^y = dr A d^ z A n. 
Then for any configuration y that satisfies the Bethe equations 

dy^S = 0, S = 1,...,K 

one derives from (^ ) 

HQG = —{drS)G , . 

HnG = -{d,„S)G ^^^^ 

i.e. G is an eigenvector of Hq and Hn,n = 1,. . . ,N with eigenvalues —drS and —dz„S, 
respectively 

Here we will need a generalization of these equations in order to get solutions with the 
asymptotic behavior of eq. (p^), that will eventually lead to the Richardson equations. 
Following [^] take ^ G t)*^ and define 

n = e<"'«>fi, 5 = S + -^lei V - ^(C, An)zn + as)ys, (43) 

from (^0|) one can derive a similar relation for the new objects 

dS AQ + drA HoQ + ^ A Hn^ = (44) 
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and then one has eigenvectors 



Foe<"'«>G = -(-^leP + 5rS)e<"'«>G 

, , 47^^ 
F,e<"'«>G = ((C,A„)-5,„5)e<"'«>G 



provided y is chosen so that 



s'^s n 



In order to obtain eigenvectors for the Hamiltonians of eqs. (35) we redefine ^ 
and take, as in the previous section, the two consecutive hmits r ioo and z,u ^ 
C ^ in G. 

After this hmiting procedure G ^ F where 

F{u,z,y) = ^ FK,Va 

K VgSk 

with 

FK,a{u,Z,y) = U/3i_i(yi,l - yi,2)vpi,2iyi,2 " ^l.s) " ' ' V(3^^ji^{yi,Ki " ^l) 

X 

X f/3jv,i {yN,l - yN,2)^f3N^2 iyN,2 - yN,z) ■ ■ ■ Vf3N,Kj^ {yN,KM - Zn) 

where ^ ^ 

vpiz) = -- - — 
z {u,P) 

and the rest of the notation is as in eq. (^) . 
One then satisfies the eigenvalue equations: 

i?„e<"'«>F=((C,A„)-a,„r)e<"'«>F 

which hold for 

T{z,y) = ^ (An, A„/)ln(z„ - z„/) - ^(A„, a^) ln(2;„ - y^) 

n<n' ri,s 

+ ^{oLs,as')\n{ys - ys') 

s<s' 

provided y is chosen so that it fulfills the Bethe equations 

dy^T+{C,as) = s = l,...,K 

or explicitly 

^i^^-$:i^^ = (e,a.) s = l,...,K. 

ys' ys n y^ 
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Note that F{u/e, z,y) is for generic z,y and u a polynomial in e thus it admits the 
expansion of eq. (^) and 

Fo{z,y) = limF{u/e,z,y) 

— e— *U 

' - an eigenvector of Hamiltonians Rn provided y is a solution of (^0|) , whose corresponding 
eigenvalue /i„ is given by (C, A„) - d^^T, i. e. 

Hn = + 2^ — 2^ — - — (51) 

In our construction we only considered the subspace of vectors in the zero total 
weight subspace, this was imposed from the very beginning, for only when restricted to 
this subspace the Hamiltonians Hn commute, defining an integrable system. However 
after the limiting procedure we arrive at Hamiltonians Rn and these have the property 
of commuting in the whole representation space. We should then produce solutions of 
arbitrary weight. 

Suppose we search for an eigenvector of the Hamiltonians Rn of total weight Aq- In 
order to get it we introduce a new insertion at zq with representation of highest weight Aq 
and vector space Vaq- Now we decompose the sum of highest weights 

L N 
^ a, = ^ A„ + Ao 

s=l n=l 

with Aq and the last L — K roots chosen to verify 

L 

s=K+l 

We reparametrize the corresponding screening variables 

{yu---^yL) = {yi,...,yK:Zo+-fi,...,zo + 'yL-K) 

and send zq to oo while keeping and 7; finite. In the limit, the old Hamiltonians Rn 
decouple from the new insertion, the Bethe equations (|50|) for s = 1, . . . ,K are unchanged 
(the new terms coming from the new insertion and screening variables vanish in the limit) 
and the eigenvectors factorize in Fq 0/ where / € Vaq has weight — Aq- As the total weight 
is zero Fq is an eigenvector of Hamiltonians Rn (provided the Bethe equations are fulfilled) 
of weight Aq. In practice we can forget about the new insertion and screening variables 



and simply solve the Bethe equations (^) and use the Bethe ansatz solution (46p7|) but 
with a new set of simple roots so that 

K N 

^a, = ^A„-Ao, (52) 

s=l n=l 

and accordingly a different number of screening variables ys- 

This ends our main results which show the deep mathematical structure underlying the 
connection of BCS-like Richardson models with Chern-Simons theory. 
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6 Conclusions 



In this paper we have shown how the BCS model, its exact solution, and integrability is 
related to the Chern-Simons theory. This has been done in two ways. First by adding a 
perturbation to the CS action on the plane by means of a one-dimensional coloured charge 
distribution, and second by taking an appropriated limit of the CS model on a torus with 
twisted boundary conditions for the WZW gauge fields. The construction is valid for any 
simple, simply connected compact Lie group G. The case of BCS corresponding to the 
choice G = SU{2). 

We find quite remarkable that the Chern-Simons theory not only gives effective de- 
scriptions of the fractional quantum Hall effect and the anyon superconductivity, but also 
of the standard BCS superconductivity. Though in the later case the formalism works, 
not in real space, but in the complex energy space. The space dimensionality enters only 
through the representations of the Lie group associated to every energy level, which for 
BCS is simply the level degeneracy. 

This work adds one more member to the series: BCS/Gaudin models — > Integrable Ver- 
tex Models — > WZW models — > Chern-Simons theory, suggesting new research directions. 
One is the possible existence of integrable models associated to higher genus Riemann 
surfaces. This is an interesting but difficult problem, see for a study of the genus 2 
case. 

Another question comes from the free field realization of the SU{2) WZW model. This 
model has two screening operators J±, where only one of them, i.e. J+, is required for 
the computation of correlators of the WZW primary fields with integer spins, while the 



second one is needed for primary fields with fractional spin 42, 43|. In the limit 



where k + 2 ^ i), the positions of the screening operators J+{u) get frozen and satisfy the 
Richardson/Gaudin equations 0]. Thus J+ can be regarded as a Cooper pair creation 
operator. The question is: what is the role, if any, of the second screening operator J_ and 
the associated primary fields with fractional spin? A conjecture is that they are associated 
to vortex operators in the BCS model since they produced additional zeros in the wave 
function. This possibility is currently under investigation. 
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Appendix 

The aim of this appendix is to give a more explicit expression of the basic equations 
(50,51) for their application to several problems in Condensed Matter and Nuclear Physics. 
To this end we shall exploit the structure of the simple Lie groups, together with their 
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representation theory |44] 



The problem that is posed is the following: given a set of complex variables {zn}n=i 
together with a set of highest weights (h.w.) {A„}^^i of a simple Lie group G and a generic 
element ^ belonging to the Cartan subalgebra f)*^ of G, find the solutions of "semiclassical 
Bethe equations" ( [50| ) for a set of K unknowns {ys}^=i, which are associated to a set of 
simple roots {as}^=i- Moreover the sets of h.w. and simple roots are related by the eq. 
(52) where Aq is the total weight of the eigenvector we look for. Once the solutions of (pO|) 



are found then eqs.(|5l|) give the eigenvalues /i„ of the "perturbed" Gaudin operators Rn 



defined in eqs. (38). 



A simple Lie Group G, with rank r, is characterized by r simple roots {ota}a=i- It 
is convenient to describe the set {as}^=i by the numbers {Na}a=i of roots of the same 
type cta- Then the label s amounts to the pair of labels {a,i) with a = 1,... ,r and 
i = 1, . . . ,Na, such that K = J2a=i ^a- Correspondingly we shall make the change of 
variables: ys = y[a,i) = ^t- With these definitions eqs.(^0|) become 

EE f^-Ef^ = (53) 

b=lj=l % n=l ^n-^i 

where a = 1, . . . , r and i = 1, . . . , Na- Of course the term i = j,a = b must be excluded in 
jSt). Similarly eqs. (|5l|) read, 



f,n = {^,An)+ }^ -1^1^ pa _^ (54) 

The various scalar products (•, •) appearing in these eqs. are well known in the theory 



of Lie algebras |44]. First of all, the scalar product of simple roots {a.a,ab) is directly 
related to the Cartan matrix Aab, which is defined as 

Aab = {aa,ab) (55) 

where is the coroot of cta given by 



< = ^ (56) 



For simply laced Lie groups (i.e. the ones in the ADE series) and in the normalization 
we have adopted, one has lotaP = 2, Va, so that roots and coroots are equal, and the 
Cartan matrix is symmetric. 

Next we have to compute the scalar products involving the h.w.'s A„. Every weight A 
is characterized by its expansion in the basis of the fundamental weight vectors {uja}a=i, 
i.e. 

r 

A = ^ A-^ u;, (57) 

a=l 
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where {A^jJ^^^^ are integers called Dynkin labels (they are non negative for h.w.v). The 
basic property of the fundamental weights is 



ab 



Using eqs(56,57,58) one can deduce 



{\OLa) 



The scalar product of two weights Ai and A2 is given by 

r 

(Al,A2) = KFabA = ^i-F-\2 

a,b=l 

where Fab is called the quadratic form matrix which is defined as 

Fab = {i^a,i^b) = {^'^)a,b^r- 



(58) 



(59) 



(60) 



(61) 



The last equality in eq. ( pi] ) follows from the relation between simple roots and fundamental 
weights, 



"a = X! ^ab^b 
6=1 

The element ^ G f)*^ will also be expanded in the basis of fundamental weights, 
Then using eqs. ( ^9|) and (|60|) we get, 



(62) 



(63) 



a=l 



a 



(e,A)=e-F.A, 

Using the formulae introduced above we can finally write the Bethe eqs.(^3|) as 

r Aft, ^ N 



(64) 



\ ^ \ ^ ^ba _ \ ' ^^n 

E^-- Ef ^ z„. - Ef 



(65) 



6=1 j=l J « n=l 

and the eigenvalues (|5^) as 



A„, ■ F- 



r Na 



a=l i=l ^" 



(66) 



where A" are the Dynkin labels of the h.w. A„. In reference [44| can be found the explicit 
expressions of the Cartan and quadratic form matrices, so the remaining data to be fixed 
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are the representations associated to the sites as weh as the weight Xq (see below for 
examples) . 

Example 1: G = SU{2) 

SU{2) is the unique simple Lie group with rank 1 (i.e. r = 1), with a simple root 
cxi = \f2 and a fundamental weight <jJ\ = ai/2 = l/\/2- The Cartan and the quadratic 
form "matrices" are given by An = 2 and Fn = 1/2. The Dynkin label of a h.w. A is 
twice the spin, i.e. = 2s, (s = 0, 1/2, 1, . . .). With the notation A^i = M (M being the 
number of Cooper pairs) and = —1/g, {g being the BCS coupling constant), eqs.(|65|) 
yield the Richardson eqs. (^) 

2 ^ 2s 1 



associated to the BCS Hamiltonian, while eqs.(p6D yield the eigenvalues of the CRS con- 
served quantities 1 11 , 12 , • 



r) Mr) 

y n'j^n ^" ^" i=l ^" 

In the BCS case, is the number of energy levels which are occupied by electrons 
pairs, Zn is twice the single particle energy (i.e. z„ = 2e„), 2s„ is the maximal number of 
pairs that can occupy that level. 

Denoting by sq the spin of the weight Aq appearing in eq. (^2|) one gets 



N 

SO = E ^« - ^ (69) 

71=1 

In the case where Sn = 1/2, Vn, the half-filling constraint (i.e. as many electrons as energy 
levels) is given by A^ = 2M, which from ( |69| ) corresponds to sq = 0. 

Before we consider more general groups, we shall discuss in more detail eq.(|52|) regarding 
the relation between the number of roots Na and the weights A„ and Aq. First of all this 
eq. can be written as 

r N 

J2 NaCXa = ^ A„ - Ao (70) 
a=l n=l 

The problem is to find the numbers Na given A^, and Aq- Actually there may 

not exists a solution of this eq. with all A^a non negative integers. 

A particular case where we know there must exist a solution is for all the A^s the 
adjoint irrep and Aq = 0. The h.w. of the adjoint representation is given by the highest 
root 6 



9 = J2 daOta (71) 



a=l 
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where da are positive integers called the marks (we use a different notation from reference. 
Ql). For SU{n) groups da = I, Va. Setting A„ = 0, Vn and Aq = in ( [70|) the solution 
is given by 

Na = daN (72) 

Let us now consider the more general case. Using eq.(B2|) in ([70[) we get 



r N 

J2N,A,a = Y.K-K (73) 

6=1 n=l 



Multiplying by the inverse of the Cartan matrix and using eq.(|6lD we finally get the desired 
expression: 

Na = T;^j:na(j:^n-X'^ (74) 
b=l \ n I 

which obviously may not yield integers values for Na- 
Example 2: G = SU{r + 1) 

The Lie group A^. = SU{r + 1) has rank r. We shall choose N insertions in the 
fundamental representation, whose Dynkin labels are given by A = (1,0, ...,0). From 
eq.(|7^) we derive then 

Xl= N -2Ni + N2 
Xl= -2N2 + N1 + NS 

(75) 

y = -2Nr + Nr-l. 



If we choose Aq = then the number of roots Na follows from eq.(74) 

JV„ = iV '-t^ (76) 

r + 1 

Hence must be a multiple of r + 1. This reflects the fact that the tensor product of r + 1 
copies of the fundamental irrep contains the identity. 

The Bethe eqs. in the case of SU{3) and arbitrary irreps is given by 

\ ^ ^ \ ^ ^^n _ a1 

E\ -E^ ^ E"^ -E^ ^ z -E^ 

A'2 r, Afi T ^ A2 ^ ' 
\ ^ ^ \ ^ ^^n _ c2 

. K2 _ ^2 E] - E^ ^,Zn- Ef ^ 

j=l J * j=l J * n=l * 

Other representations and groups can be worked out similarly. 
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